Disappearance of Quantum Chaos in Coupled Chaotic Quantum Dots 
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Statistical properties of the single electron levels confined in the semiconductor (InAs/GaAs, 
Si/Si02) double quantum dots (DQDs) are considered. We demonstrate that in the electronically 
coupled chaotic quantum dots the chaos with its level repulsion disappears and the nearest neighbor 
level statistics becomes Poissonian. This result is discussed in the light of the recently predicted 
"huge conductance peak" by R.S. Whitney at al. (Phys. Rev. Lett. 102, 186802 (2009)) in the 
mirror symmetric DQDs. 
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The problem of chaos in nanostructures has a relatively 
long history just since these structures effectively entered 
science and technology. The importance of this problem 
is connected with wide spectrum of the transport phe- 
nomena. The problem was actively studied in the last 
two decades [T]-[4] . One of the main results of these stud- 
ies, based on the classical and semi-classical approaches, 
is that the various transport phenomena sensitively de- 
pend on the geometry of the objects and, first of all, on 
the symmetry: Right - Left (RL) mirror symmetry, up- 
down symmetry and preserving the loop orientation in- 
version symmetry important in the presence of the mag- 
netic field [5 ,E] (see also earlier paper [7 ). 

There is another, actively studied in various fields of 
physics, aspect in some sense complimentary: Quantum 
Chaos with its inalienable quantum character including, 
first of all, well known Nearest Neighbor level Statistics 
(NNS) which is one of standard quantum-chaotic test. In 
our recent paper [8l , we investigated the NNS for various 
shape of the single quantum dots (SQD) in the regime 
of the weak confinement when the number of the levels 
allows to use quite sufficient statistics. Referring for de- 
tails to [8], we briefly sum up the main conclusions of 
[8]: SQDs with RL or up-down mirror symmetries have 
a Poisson type NNS whereas a violation of both of these 
symmetries leads to the Quantum Chaos type NNS. 

In the present paper we study a quantum chaotic prop- 
erties of the double QD (DQD). By QD here we mean 
the three dimensionally (3D) confined quantum object, 
as well its 2D analogue - quantum well (QW). In three 
dimensional case we effectively use, as in previous paper 
Ref. [8], an assumption of the rotational symmetry of 
QD shape. We follow the approach developed in Ref. [8 
for the calculation of the levels of the QD. QDs embedded 
in substrate are considered in the single sub-band effec- 
tive mass approximation with envelope functions supple- 
mented by the energy dependence of the electron effective 
mass and non-parabolic effect. This approach is in good 
agreement with the experimental data and previous cal- 
culations in the strong confinement regime 0[9]. Here, in 
the regime of weak confinement, as in Ref. [8 , we also do 
not consider Coulomb interaction between electron and 



hole: Coulomb effects are weak when the barrier between 
dots is thin leading to the strong interdot tunneling and 
dot sizes are large enough. In these circumstances, stud- 
ied in detail in Ref. [10 (see also for short review a mono- 
graph [11 ), one may justify disregard of the Coulomb 
effects. The physical effect, we are looking for, has place 
just for thin barriers; to have sufficient level statistics, we 
need large enough QDs (> 100 nm for InAs/GaAs QW). 
Thus, we consider tunnel coupled two QDs with substrate 
between, which serves as barrier with electronic proper- 
ties distinct from QD. Boundary conditions for the single 
electron Schrodinger equation for wave function ^(r) and 
its derivative l/m*(n, V)^(r) on interface of QD and the 
substrate are standard. We take into account the mass 
asymmetry inside as well outside of QDs [8]. To avoid 
the complications connected with spin-orbit coupling, s- 
levels of electron are only considered in the following. We 
would like to remind that the selection of levels with the 
same quantum numbers is requisite for study of NNS and 
other types of level statistics. 

Whereas at the large distances between dots each dot 
is independent and electron levels are twofold degenerate, 
expressing the fact that electron can be find either in one 
or in the other isolated dot, at the smaller interdot dis- 
tances the single electron wave function begins to delocal- 
ize and extends to the whole DQD system. Each twofold 
degenerated level of the SQD splits by two, difference of 
energies is determined by the overlap, shift and transfer 
integrals [12 . Actually, due to the electron spin, there 
is fourfold degeneracy, however that does not change our 
results and below we consider electron as spinless. Note 
that the distance of removing degeneracy is different for 
different electron levels. This distance is larger for lev- 
els with higher energy measured relative to the bottom 
quantum well (see Fig. [I] below). By the proper choice 
of materials of dots and substrate one can amplify the 
"penetration" effects of the wave function. 

The single electron energy levels are notated by Ei, 
i = 0, 1, . . . , A^. Energy differences between neighboring 
levels are AEi = Ei — i = 1, 2, . . . , A^. The normal- 

ized distribution function R{AE) is calculated following 
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FIG. 1. (The left figure) The band gap model for InAs/GaAs 
heterostructures. (The right figure) The band gap model for 
Si/Si02 heterostructures. The strength value of band gap 
and strain induced (0.21 eV for InAs/GaAs, negligible for 
Si/Si02) potentials and effective masses are shown |8|l9j. CB 
(VB) means conductive band (valence band). 



the procedure of Ref. [8] : 

R{AEj) = Rj=Nj/Hae/N, 



j = l,2,...,M, 



where N = Nj represents total number of considered 
levels, Hae = {^En — AEi)/M is the energy interval 
which we obtained by dividing the total region of energy 
differences by M bins. Nj is the number of energy dif- 
ferences which are located in the j-th bin. 

We find the distribution functions R{AE) using the 
smoothing spline method. If Ri, i = 1, . . . , M, are cal- 
culated values of the distribution functions correspond- 
ing to AEi^ the smoothing spline is constructed by giv- 
ing the minimum of the form ^f^i{Ri — R{AEi))'^ + 
/ R"{AEYd{AE)/\. The parameter A > is control- 
ling the concurrence between fidelity to the data and 
roughness of the function sought for. For A ^ oo one 
obtains an interpolating spline. For A ^ one has a lin- 
ear least squares approximation. Below we display some 
of our results for semiconductor DQDs. In Fig. [l] we 
present schematically the band gap model (see Refs. [8l|9] 
for details) with used in this paper materials and corre- 
sponding parameters of QDs. Fig. |2] shows distribution 
function for two Si/Si02 QDs of the shape of the 3D el- 
lipsoids with a cut below the major axis. Isolated QD of 
this shape is strongly chaotic. It means that distribution 
function of this QD can be well fitted by Brody formula 
[13] with the parameter which is close to unity [8^. We 
see that the corresponding up-down mirror symmetric 
DQD shows Poisson-like NNS. Note that these statistics 
data involved 300 confined electron levels, which filled 
the quantum well from bottom to upper edges. We con- 
sidered the electron levels with the orbital momentum 
/=0, as was mentioned above. The orbital momentum 
of electron can be defined due to rotational symmetry of 
the QD shape. 

In Fig. |3] SQD (2D quantum well) without both type 
of symmetry reveals level repulsion, two tunnel coupled 
dots show the level attraction. From the mirror sym- 
metry point of view, the chaotic character of such single 
object is due to the lack of the R-L and up-down mirror 
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FIG. 2. (The upper figure) The cross section of DQD shape 
is shown (sizes are given in nm). The electron wave func- 
tion of the ground state is shown by the contour plot. (The 
lower figure) Distribution functions for energy differences of 
the electron neighboring levels in Si/Si02 single QD (solid 
line) and DQD (dashed line). The coefficient of the spline 
smoothing is equal 6. 



symmetries. The symmetry requirements in this case, for 
the coupled dots are less restrictive: presence of one of 
the mirror symmetry types is sufficient for the absence 
of quantum chaos. 

Dependence NNS on the interdot distance shows a 
gradual transition to the regular behavior with interme- 
diate situation when Poisson-like behavior coexists with 
chaotic one: they combine but the level attraction is not 
precisely Poisson-like. Further decreasing distance re- 
stores usual Poisson character (see Fig. |3|. 

Fig. [4] shows how the degeneracy gradually disappears 
with the distance h between QDs in InAs/GaAs DQD. 

Finally, we would like to show the disappearance of 
the Quantum Chaos when chaotic QW is involved in the 
"butterfly double dot" [6] giving huge conductance peak. 
Fig. [5] shows the NNS for chaotic single QW of [6 by 
dashed line. Mirror (up-down and L-R) symmetry is vi- 
olated. The NNS for a L-R mirror symmetric DQW is 
displayed by solid line in Fig. |5j It is clear that Quantum 
Chaos disappears. 
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FIG. 3. (The upper figure) The DQW shape is shown (sizes 
are in nm). (The lower figure) Distribution functions for en- 
ergy differences of the electron neighboring levels in the 2D 
InAs/GaAs DQW calculated for various distances b between 
QWs. The coefficient of the spline smoothing is equal 4. Dis- 
tribution functions of single QW is shown. 



We conjecture that the above mentioned peak in con- 
ductance of Ref. [6] and observed here a disappearance 
of Quantum Chaos in the same array are the expres- 
sion of the two faces of the Quantum Mechanics with its 
semiclassics and genuine quantum problem of the energy 
levels of the confined objects, despite the different scales 
(what seems quite natural) in these two phenomena (sev- 
eral micrometers and 10 - 100 nm, wide barrier in the 
first case and narrow one in the second). We have to 
emphasize here that the transport properties are mainly 
the problem of the wave function whereas the NNS is 
mainly the problem of eigenvalues. Similar phenomena 
are expected for the several properly arranged coupled 
multiple QDs and QD superlattices. In the last case, 
having in mind, for simplicity, a linear array, arranging 
the tunnel coupling between QDs strong enough, we will 
have wide mini-bands containing sufficient amount of en- 
ergy levels and the gap between successive minibands will 
be narrow. Since the levels in the different minibands are 
uncorrelated, the overall NNS will be Poissonian indepen- 
dently of the chaotic properties of single QD. We would 
like to remark also that our results have place for 3D as 
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FIG. 4. (The upper figure) Doublet splitting of single 
electron levels dependence on the distance b between QDs in 
InAs/GaAs DQD. The ground state {E = 0.23 eV) level split- 
ting is expressed by dashed line. The solid line corresponds to 
doublet splitting of a level which is close to upper edge of the 
quantum well {E = 0.56 eV). The shape of DQD is the same 
as in Fig. [5] (The lower figure) The electron wave functions 
of the doublet state: the ground state (left) and first excited 
state (right), are shown. 



well as for 2D quantum objects. It is important to notice 
that the effect of reduction of the chaos in a system of 
DQD could appear for interdot distances larger than con- 
sidered, for instance in figures [3j if an external electrical 
field is applied. By properly designed bias, the electric 
field will amplify wave function "penetration", effectively 
reducing a barrier between QDs. 

In conclusion, we have shown that the tunnel coupled 
chaotic QDs in the mirror symmetric arrangement have 
no quantum chaotic properties, NNS shows energy level 
attraction as should to be for regular, non-chaotic sys- 
tems. These results are confronted with the huge con- 
ductance peak found by the semiclassical method in Ref. 
[6 . We think that our results have more general appli- 
cability for other confined quantum objects, not only for 
the quantum nanostructures, and may be technologically 
interesting 

We thank Igor Bondarev for useful discussions. 
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FIG. 5. Distribution functions for electron neighboring levels 
in InAs/GaAs single QW (dashed line) and DQW (solid hne). 
The coefficient of the spline smoothing is equal 4. Shape of 
DQW is shown in the inset. The electron wave function of 
the ground state is shown by the contour plot in the inset. 
Data of the statistics include 200 first electron levels. 
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